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Abstract 

A canonical quantization scheme is represented for a quantum sys- 
tem interacting with a nonlinear absorbing environment. The envi- 
ronment is taken anisotropic and the main system is coupled to its 
environment through some coupling tensors of various ranks. The 
nonlinear response equation of the environment against the motion 
of the main system is obtained. The nonlinear Langevin-schrodinger 
equation is concluded as the macroscopic equation of motion of the 
dissipative system. The effect of nonlinearity of the environment is 
investigated on the spontaneous emission of an initially excited two 
level-atom imbedded in such an environmrnt. 
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1 Introduction 



The simplest way of describing a damped system in classical dynamics is 
by adding a resisting force, generally velocity-dependent, to the equation of 
motion of the system. Frequently the magnitude of the resisting force may 
be closely presented, over a limited range of velocity, by the law fd = av n , 
where v is the velocity of the damped system and a and n are constants. 
For example for the friction force n — 0, viscous force n = 1 and for high 
speed motion n = 2 [TJ. Such an approach is no longer possible in quantum 
mechanics, because one can not find a unitary time evolution operator for 
both the states and the observables, consistently. 

In order to take into account the dissipation in a quantum system, there are 
usually two approaches. The first approach is a phenomenological way, by 
which the effect of dissipation is taken into account by constructing a suitable 
Lagrangian or Hamiltonian for the system [2j [3] . Following this method the 
first Hamiltonian was proposed by Caldirola jl] and Kanai [5] and afterward 
by others [21 E] • There are difficulties about the quantum mechanical solu- 
tions of the Caldirola-Kanai Hamiltonian. For example quantization using 
this way violates the uncertainty relations or canonical commutation rules. 
The uncertainty relations vanishes as time tends to infinity [8]- jllj . 

The second approach is based on the assumption that the damping forces 
is caused by an irreversible transfer of energy from the system to a reservoir 
[T2[ [13]. In this method , modeling the absorptive environment by a col- 
lection of harmonic oscillators and choosing a suitable interaction between 
the system and the oscillators, a consistent quantization is achieved for both 
the main system and the environment [14]- 126|. In the Heisenberg picture, 
one can obtain the linear Langevin-schrodinger equation, as the macroscopic 
equation of motion of the main system. [T4l [T5] . 

In the present work, following the second approach, a fully canonical 
quantization is introduced for a system moving in an anisotropic non-linear 
absorbing environment. The dissipative system is the prototype of some 
important problems which the present approach can be applied to cover such 
problems straightforwardly. 

The paper is organized as follows: In section 2, a Lagrangian for the total 
system (the main system and the environment) is proposed and a classical 
treatment of the dissipative system is achieved. In section 3, the Lagrangian 
introduced in the section 2 is used for a canonical quantization of both the 
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main system and the non-linear environment. In section 4, the present quan- 
tization is used to investigate the effect of the nonlinearity of the environment 
on the spontaneous decay rate of an initially excited two-level atom embedded 
in the absorbing environment. Finally, the paper is closed with a summary 
and some concluding remarks in section 5. 

2 Three-dimensional quantum dissipative sys- 
tems 

When an absorbing environment responds non-linearly against the motion of 
a system, the non-linear Langevin-Schrodinger equation is usually appeared 
as the macroscopic equation of motion of the system. As an example, when 
the electromagnetic field is propagated in an absorbing non-linear polarizable 
medium, the vector potential satisfies the non-linear Langevin-Schrodinger 
equation. In this section, the motion of a three-dimensional system in the 
presence of an anisotropic non-linear absorbing environment is classically 
treated. For this purpose , the environment is modeled by a continium of 
three dimensional harmonic oscillators labeled by a continuous parameter u. 
The total Lagrangian is proposed as 



which is the sum of three pars. The part L e is the Lagrangian of the envi- 
ronment 



where X(a>, t) is the dynamical variable of the oscillator labeled by uj. The 
second part L s in ([1]) is the Lagrangian of the main system. Taking the 
system as a particle with mass, m, moving under an external potential V(q), 
one can write 



The last part Li nt in the total Lagrangian (pQ) is the interaction term between 
the system and its absorbing environment and includes both the linear and 



L(t) — L e + L s + Li n t. 



(1) 
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nonlinear contributions as follows 

Lm = J™ duf!p(u))?(t)X3{u,t) 

p OO P OO 

+ dcuJ o du'fgl(u,u')q\t)X^u,t)X k (uj',t) 

/•oo r oo /'oo 

+ / dcj / dw' / du;" #?2,(w, c/, £*/')?*(*) X J '(w, £)X fe (u/ 3 t)X l (u", t) + . . 

JO JO •/ 

where f^ 2 \ f^ 3 \ ■ ■ ■ are the coupling tensors of the main system and 
its environment. As it is seen from (Tj0) the coupling tensor describes 
the linear contribution of the interaction part and the sequence f( 2 \ f^ 3 \ .... 
describe, respectively, the first order of the non-linear interaction part, the 
second order of the non-linear interaction part and so on. The interaction La- 
grangian (j4]) is the generalization of the Lagrangian that previously has been 
applied to quantize the electromagnetic field in the presence of anisotropic 
linear magnetodielectric media |27j . 

The coupling tensors f^\f^ 2 \f^ 3 \ ■■■ in (j3J) are the key parameters of this 
quantization scheme. As it will be seen, in the next section, the susceptibil- 
ity tensors of the environment (of the various ranks) are expressed in terms 
of the coupling tensors. Also the noise forces are obtained in terms of the 
coupling tensors and the dynamical variables of the environment at t = —oo. 

2.1 The classical Lagrangian equations 

the classical equations of motion of the total system can be obtained using 
the principle of the Hamilton's least action, 5 J dt L(t) = 0. These equations 
are the Euler-Lagrange equations. For the dynamical variables, X(u>,£), the 



Euler-Lagrange equations are as 



d 



5L 



5L 



1,2,3 



X i (u,t) + tSX i (u,t) = ?{t)f£\u) 



+ 



+ 



(2), 



dJ I dco" qi(t) f%(u, co', co") + f^', co, co") 



(3) 



+ f%(co',co",u;) X k (co',t)X(u",t) + 



(5) 



Also the Lagrange equations for the freedom degrees of the main system are 
obtained as follows 



d 



6L 



5L 



where 
Ri(t) 

+ 



dt \5(q t (t))J %(t)) 
=> mq(t) + VV(q) = -R(t) 



dufg\u)Xi(u,t) + I du I dco'm (u,u')X^u,t) X k (u',t) 







1,2,3 



(6) 



dco 



dw'fl%{uj,uj\Lo")X^u,t) X k (co,t) X l (co",t) + - 



(7) 



In Eq. ([6]) — R(t) is the force exerted on the main system due to its motion 
inside the absorbing environment. It will be seen that the force — R(t) can be 
separated into two parts. One part is the damping force which is dependent 
on the various powers of the velocity of the main system. The second part is 
the noise forces which has sinusodial time dependence. Both the damping and 
the noise forces are necessary for a consistent quantization of a dissipative 
system. Without the noise forces the quantization of a dissipative system 
encounter inconsistency. According to the fluctuation- dissipation theorem 
the absence of any of these two parts leads to the vanishing of the other part. 
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3 Canonical quantization 



In order to represent a canonical quantization, the canonical conjugate mo- 
menta corresponding to the dynamical variables X(o;, t) and q should be 
computed using the Lagrangian (JTJ). These momenta are as follows 

SL 

Qi(u,t)= =Xi{u,t) i = 1,2,3 (8) 

6(Xi(u,t)) 

S L 

Pi{t) = ^-y = mft + Ri{t) i = l,2, 3. (9) 

Having the canonical momenta, both the dissipative system and the envi- 
ronment can be quantized in a standard fashion by imposing the following 
equal-time commutation rules 

, Pj (t)] = iHSij (10) 

[Xi(u, t) , Qjiu', t)} = ihS^u - u') (11) 

Using the Lagrangian ([1]) and the expressions for the canonical momenta 
given by (jSJ) and ([HD, the Hamiltonian of the total system clearly can be 
written as 

H(t) = [p(t) ~ R(t)]2 + V(q) + I [°° du \Q 2 (uj, t) + o; 2 X 2 (u;, t)l(12) 
2m 2 Jo 1 J 

where the cartesian components of R(£) is defined by ((7j). The Hamiltonian 
(|12p is the counterpart of the Hamiltonian of the quantized electromagnetic 
field in the presence of magnetodielectric media Using the commu- 
tation relations (jTUl) . (II ip and applying the total Hamiltonian (fT2|) . it can 
be shown that the combination of the Heisenberg equations of motion of the 
canonical variables X(a;,£) and Q(u,t) leads to the Eq.Q. Similarly, one 
can obtain Eq. (JH]) as the equation of motion of q(t) in the Heisenberg picture. 
Let us introduce the annihilation and creation operators of the environment 
as follows 



6i(w, t ) = \ji^- l^M^, t) + iQi(u, t)] . (13) 

From the commutation relations (|lip it is clear that the ladder operators 
bi(uj,t) and b\(u,t) obey the commutation relations 

bi(u,t),btiu',t)] = 5 i3 5{uj - uj ! ) (14) 
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The Hamiltonian (1121) can be rewritten in terms of the creation and annihi- 
lation operators bi(u,t) and b\(u,t) as follows 



H 



(p - mr 

2m 



+ V(q) + H v 



where 



3 . 

H m = Y^ / duhuj b\(uj,t)bi(u,t) 

i=l 



(15) 



(16) 



is the Hamiltonian of the absorbing environment in the normal ordering form 
and 



Ri(t) 



/ fl 



+ duo I duj f 



h 



= uj') [&>, t)b k (uo', t)+ + b](u, t)bl(u', t) 



2fuw 

h- 
j 



+bj(u, t)b\{u', t) + bUto, t)b k {uj', t) 



+ 



(17) 



are the cartesian components of the operator R, where the summation should 
be done over the repeated indices. 



3.1 The response equation of the environment 

The response equation of the absorbing environment is the base of separating 
the force — R(t) , in the right hand of ([6]), into two parts, that is, the damping 
force and the noise force. If Eq.(jSJ) is solved for X(o>, t) and then, the obtained 
solution is substituted into the definition of R(t) given by (j7|), one can obtain 
the response equation of the environment. The differential equations ([3]) are 
a continuous collection of coupled non-linear differential equations for the 
dynamical variables X(w, t). The exact solution of this equation is impossible 
unless an iteration method to be used. For simplicity here we apply the 
first order of approximation and neglect the terms containing the coupling 
tensors f( 2 \ f^ 3 \ ... in the right hand of (j^D and write the solution of Eq. (!5]) , 
approximately, as 

X(u;,t) = X N (u;,t)+ f dt'^^- (/^(u) • q(t'), 

(18) 
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where (f^)jj(u>) = fjj\u) and Xjv(u;,£) is the solution of homogeneous 
equation Xat(co>, t) + u 2 Xn{u, t) = 0. In fact Xtv(^, t) is asymptotic form of 
X.(oj,t) for very large negative times and can be written as 

X m (co,t) = [br(u)e- wt + bf(u)e^} (19) 

where b rn (u) and b^ n (u) are some time independent annihilation and creation 
operators which obviously satisfy the same commutation relations (1141) . The 
approximated solution ffl8|) yields the response equation of the environment, 
such that, the susceptibility tensors appearing in it, satisfy the various sym- 
metry properties reported by the literature [30] . 



Now substituting X(a;,i) from ([T81 in ([7j), the response equation of the non- 
linear absorbing environment is found as follows 

R(t) = R« + R^ + .... 

J — oo 

^ (2) (t) = f°° eft' / + °° <ft"xg>(i - t\t - t'%(t') Qk(t") + Rpi{t) 

J —oo J —oo 

(20) 

where x 1S the susceptibility tensor of the environment in the linear regime 
and is defined by 

X { ${t) = \ J ° Ul (21) 

[ t < 

(2) 

and Xi/k causes the first order of the nonlinearity of the response equation, 
where for for ti,t 2 > are given by 

( 2 )(+ +\ f°° a f°° a sinuJ ih sinw 2 t a {2) {1) (1) 

J JO JO U3\ U2 

(22) 

and xlfk * s zero f° r t 1 ,t 2 < 0. In (121"]) and fl22l the summation should be 
done over the repeated indices m, n. From the definition (12T1) it is clear 
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that is a symmetric tensor, Xij = X% ■ There are also some symmetry 
features for the non-linear susceptibility tensors of the various orders. These 
symmetry properties can be satisfied by imposing some conditions on the 
coupling tensors f( 2 \ f( 3 \ ... . For example the susceptibility tensor x ■> 
should satisfy the symmetry property [30] 



xg!(«i,*a) = xS(Mi) (23) 

where is fulfilled provided that the coupling tensor obey the symmetry 
condition 

f$M = f!$V,<<>), (24) 
Similarly inserting the approximated solution X(c<j, t) from (fl8|) into ([7]), one 
can obtain the (n — l)'th susceptibility tensor of the environment in the 
non- linear regime for ti, t 2 , ...t n > as the following 

(n) /. . 4 x f°° j f°° j f°° j sin^iti smuj 2 t 2 smu n t n 
Xi i 1 ...i n {ti,t 2 ,...,t n ) = / dwi / du 2 .... duj n .... 

JO JO JO Ll>i lo 2 u n 

x fin j2 .... jn (^i^2, - w„) /iS^Wi) fQ 2 (uJ 2 ).... fQ n M (25) 

and Xi^...i„(^i) ^2, ^n) is identically zero for ti, £ 2 , ...t n < 0. the Susceptibil- 
ity tensor Xui...i n (^i)^2, ■■■,t n ) should satisfy the symmetry relations [30] 

Xi ^2) ■••) ^fcj •••£/> •■•) tn) = Xi i 1 ...ii,...i k ,...,i n {^li ^2; •■•) tli ■■■tki ■■■■> ^n) 

(26) 



where this symmetry relation is clearly fulfilled by imposing the symmetry 
conditions 

fihh J fc ,...,ii,..,inK'^2, ...,w fc , w n ) 

= 4hk ^...,^,...^(^1^2,-,^,...,^,..., u n ) (27) 

on the n'th coupling tensor in the interaction Lagrangian (J3J). 

In Eq. d20D R$(*) and R^(t) are the noise forces in the linear regime 
and the first order of non-linearity , respectively, and using the symmetry 



9 



relation (|24p are obtained as 



ij 

oo 



p OO P oo 

+ r du r dj flu*,, j) &2 {j) f df i in u ' {t - t>] 



X 



UJ' 

X n N {uJ,t)q3{t') + q\t')X™{uj\t)\ (28) 



where the summation should be done over the repeated indices and Xjv(w, t) 
is the asymptotic solution (|T9j) . 

It is remarkable that for some known susceptibility tensors X \ X ■> •••■> X > 
the coupling tensors / , / , -.-.f^ satisfying the definitions (j2B . (!22l) and 
(1251) are not unique. In fact if the coupling tensors /W, f( 2 \ f( n > sat- 
isfy (l2Tl ) . (l22l) and ( 1251) for the given susceptibility tensors, also the coupling 
tensors /'« , /'< 2 ) /'<") defined by 

/■'(i) _ f (i) A . 

f'(") - f(") 4. . 4. . 4. . roq^i 

Ji iiia. ...t n J* hh-jn - n ~n3i^2j2----- n -in]n \^) 

satisfy ( l2~Tj) .( l22"l) and ( |25i) . where A is an orthogonal matrix A im A m j = 5ij. 
The various choices of the coupling tensors f( 2 \ which is related 

to each other by the orthogonal transformation (1291) do not change the physi- 
cal observables. The commutation relations between the dynamical variables 
of the total system remain unchanged under the orthogonal transformation 
(T2T?j) . For example in the next section it is shown that the decay rate of an 
initially excited two-level atom, embedded in a non-linear absorbing environ- 
ment, are independent of the various choices of the coupling tensors which is 
related to each other by the transformation (1291) . 

Now combination of the response equation fl20l and equation ([6]) yields the 
non-linear Lagevin-schrodinger equation 



mq l (t) + / dt'x$ ) (t-t')qi(t') + 

J —CO 

dt' / dt" X ^l(t-t',t-t"W(t')q k (t") + .... 

-oo J —oo 

dV ® -flg> t ( t )-J# iW+ ... (30) 



dq l 
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as the macroscopic equation of motion of the main system in the anisotropic 
non- linear absorbing environment. The velocity dependent terms in the left 
hand of this equation are the damping forces exerted on the main system. 
The noise forces — R^i, —R^\, ... in the right hand of fl30l) are necessary for 
a consistent quantization of the dissipative system. As a realization, if this 
quantization method would be applied for the electromagnetic field in the 
presence of an absorbing non-linear dielectric medium, the vector potential 
would satisfy the equation (|30|) . In that case, the tensors X** 2 \ •••• would 
play the role of the electric susceptibility tensors and — Ryv\ — ••• would 
be the noise polarization densities of various orders. 



4 The effect of nonlinearity of the environ- 
ment on the spontaneous emission of a two- 
level atom imbedded in an absorbing envi- 
ronment 

In this section the effect of non-linearity of the absorbing environment is 
investigated on the spontaneous emission of a two-level atom embedded in 
such an environment. To calculate the spontaneous decay rate of an initially 
excited two-level atom, the quantization scheme in the preceding section is 
used and the theory of damping based on the density operator method is 
applied [31]. Neglecting the second power of the operator R in (Tl5|) the 
Hamiltonian of the total system can be written as 

H=H + H' 

2 3 

H = H S + H m = ^- + V{q) +J2 [duhw bl(uo)b t (u) 

H' = — p ■ R (31) 

Let us suppose the main system is a one electron atom with two eigenstates 
|1) and |2) correspond to two eigenvalues E\ and E 2 , respectively (E 2 > Ei). 
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the Hamiltonian (1311) can now be rewritten as [3T] , |32j . 



H = H + H' 

3 r 

i=i 

H' = % mcu R • [da — d*a^ 



E 2 — Ei 

h 



(32) 



where m is the electron mass of the atom, a = |1)(2| ,a^ = |2)(1| , are 
the Pauli operators and d = (l|r|2), where r is the position vector of the 
electron with respect to the center of mass of the atom. Dropping the energy 
noncoserving terms correspond to rotating wave approximation and regarding 
the relation (fl7|) . the interaction term H' up to the first order of nonlinearity 
in the interaction picture is expressed as 



m) = 



iHnt 

e^- H'{0) e 



du 



di<r b){u) e l( ™ }t - H.C 



+imu> / duj 



dJ 



h 



uuo' L 



2v^ 



+d l ab 3 (Lo)bl(to')e~ l{uJ0+u " UJ ' )t + diabl^bk^e-^ '^^ 1 
-<%a%(u)b k (u/)e< u °- u - u ' )t - d^<r%(u)bl(u/)e< u °- u+w ' )t 
-d^^bU^hiu')^^ ^-^ 



(33) 



where the symmetry relation (|23|) has been used. Let the combined density 
operator of the atom together with the environment is denoted by psr in the 
interaction picture. Then, the reduced density operator of the atom alone, 
denoted by ps, is obtained by taking the trace of psr with respect to the 
coordinates of the environment, that is ps = Ttr[psr\- Since it is assumed 
that H'j(t) is sufficiently small, according to the density operator approach for 
the damping theory (31], the time evolution of the reduced density operator 
ps is the solution of equation 

Ps(t) = ~Tr R [H'j(t) , p s (0) ® p R (0)] 

~Tt r J* , [H'j(t') , p s (t)®p R (0)}} (34) 
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up to order of Hf, where Pr(0) is the density operator of the environment 
at t = 0. In this formalism the environment is taken in equilibrium. Also 
the Markovian approximation has been applied replacing ps{t') by ps(t) in 
the integrand in Eq. fl34l) . 

To calculate the spontaneous emission of the atom, the initial states of the 
atom and the environment are taken as 



P*(0) = |0)<0| 



p s (0) = |2)(2| 



(35) 



where |0) is the vacuum state of the environment. Now substituting H'j(t) 
from fl3"31 into (1511) and regarding f[3"5l) the time evolution of the reduced 
density operator ps is obtained as 



Ps 



A JO U! 



m 2 UQ f°° du 



2h Jo uj 



h 



— d* fg\uj) /£» d l op s {t)J I dt' cos(a; - uj )(t - t') 

id J JO 



r(l), 



2, ,2 



m uj, 



dwi I du 2 ^— d h f-jj^uh, uh) f^j 2j2 (u 2 ,u 2 ) d i2 



(2) 



X 



t dt' apsit) ae "^.(i+i') + H C 
Jo 



duji 



oo I 

duio 



d *h /Sxii(^i^i) 4^2(^2,^2) di 



(2) 



UJ\UJ 2 
rt 



x [o*Ps(t) a + cjps{t) a" 1 "] / dt' cos uj (t — t') 



1 2 2 



00 "j 



(2) 



xaps{t) j dt' cos(u>o — uj\ — u 2 )(t — t') 

1 „ p (2) 



2, ,2 



X 



/ df e -^(*-*') <To*p s (t)+H.C 
Jo 
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m 2 ojl f°° 



/ duji / 

JO JO 



du 2 -^— d* h MUuJuUJi) f£] 2j2 {u 2 ,uj 2 ) d 



u)\U) 2 



■12 



I* dt' e — Ps (t)a ] a + H.C 

JO 



/o 

2, ,2 



m uj. 



roo r-oo I , 2 > 

Jo ^ io d * h ^uU^ 1 '^) fUihi^^V d i2 



UJIUJ2 



X 



JO 



where the the repeated indices implies that the summation should be done 
over them. Then, the equation of motion of the matrix elements psn = 
(l|ps|l) , P522 = (2|ps|2) and p S n = p* S2 i = (MPsl' 2 ) now is obtained as 



Ft Jo 00 jo 



ft 

+m 2 w 2 / / eL; 2 d* h f- l lj2 (uJi,uj 2 ) /y iia (wi,w 2 ) 4 

JO JO U)\UJ 2 
ft 



^PSTl(f) [ dt' cos( u - ^ - UJ 2 )(t - t') 
JO 



£ j( ^JO ^ 2_ ~~ ^ fifl^.^l) ./'S ,,(•• -"\ ^ /^22(0 



x[P522(*) -Psn(t)] / d* 7 coso;o(*-0 

JO 



(36) 



m 2 o;n Z" 00 rfo; 



PS22 = 



h 



°~ r — d* f^iuj) flf{u) d l p S22 (t) I dt' cos(u - UJ )(t - 
JO U JO 



-m 2 uj, 



/ du x 
Jo Jo 



du 2 -^— d* f^ ) hi2 {u 1 ,u 2 ) fi!L 2 (uJi,uJ 2 ) d 



001UJ2 



■1 2 



xps22(t) f dt' cos(u - U) X - UJ 2 )(t - t') 

JO 



duji I du 2 -^- d*^ /g^K,^) ti 2 > 2j2 (uj2,uj 2 ) d i2 p S22 (t) 

JO LO\U) 2 

t 



c(2) 



x[PS22(t) - psu(t)] [ dt' COSUJ (t-t') 
JO 



(37) 
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z Jo a; JJ 



TV? Up 1 f°° f°° I 

xp 522 (t) fdt'e-^^ (38) 



For sufficiently large times the integrals appeared in the equations ( 1361) and 
( |37|) can be approximated by 

1 ft 

— / dt' COs(uJ — UJ )(t — t') ~ 5(u — Ll>q) 



■K JO 

1 f l 



7T JO 



— / dt' cos(u;o — ijOx — ^2}{t — t') ~ 6(lu — ui\ — u 2 ) 



1 r f 

- / dt' cos <j (t - *') ~ <$M)) = (39) 

7T JO 

Hence the time evolution of the matrix elements psn and ps22 for sufficiently 
large times is reduced to 

PSll = T p S 22 PS22 = ~r P522 (40) 

where 



,2, ,2 ' 1 



+vrm u; ^ do; _ — , d* /^(w, w - w„) w - u ) d h 

(41) 

is the decay rate of the spontaneous emission of the initially excited two 
level atom up to the first order of nonlinearity. The first term in (l4"Ii) is the 
decay rate in the absence of nonlinearity effects and the second term is the 
first contribution related to nonlinear effects of the environment. It may be 
noted from (HUT) that Psn + PS22 = which implies the conservation of the 
probability. An important point is that the decay rate T is invariant under the 
various coupling tensors which is related to each other by the transformation 
). This should be so, because the decay rate T is a physical observable. 
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5 Summary 

A fully canonical quantization of a quantum system moving in an anisotropic 
non-linear absorbing environment was introduced. The main dissipative sys- 
tem was coupled with the environment through some coupling tensors of 
various ranks. The coupling tensors have an important role in this theory. 
Based on a response equation, the forces against the motion of the main 
system were resolved into two parts, the damping forces and the noise forces. 
The response equation of the environment was obtained using the Heisenberg 
equations describing the time evolution of the coordinates of the system and 
the environment. Some susceptibility tensors of various ranks were attributed 
to the environment. The susceptibility tensors in the linear and non-linear 
regimes were defined in terms of the coupling tensors of the system and its 
environment. It was shown that, by imposing some symmetry conditions on 
the coupling tensors, the susceptibility tensors obey the symmetry proper- 
ties reported in the literature. A realization of this quantization method is 
the quantized electromagnetic field in the presence of a non-linear absorbing 
dielectric. Finally the effect of the nonlinearity of the environment was in- 
vestigated on the spontaneous decay rate of a two-level atom imbedded in 
the non- linear environment. 
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